Cryptoferromagnetism in Superconductors with Broken Time-Reversal-Symmetry 
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The cryptoferromagnetic state (the state with intrinsic domain structure) in superconducting fer- 
romagnets subjected to external magnetic field is studied theoretically. Ferromagnetism originates 
either from electron spin or the intrinsic angular momentum of Cooper pairs (chiral p-wave super- 
conductors like S^RuCm). The phase transitions towards the Meissner and the mixed states are 
investigated, and the magnetic phase diagrams are obtained. Cryptoferromagnetism, as a form co- 
existence of superconductivity and ferromagnetism, can be detected by observation of magnetization 
curves predicted in the present analysis. 

PACS numbers: 74.25.Ha, 74.90.+n, 75.60.-d 
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In recent years numerous experimental evidences of 
superconductivity-ferromagnetism coexistence in various 
materials were reported 0, H, 0, 0, 0, 0]. Two types 
of such coexistence are possible: (i) The phase transi- 
tions to the ferromagnetic and the superconducting (SC) 
states occurs at different temperatures, so the coexistence 
starts below the lower from the two transitions. Ruthcno- 
cuprates [l[ belong to this type: the superconductivity 
onset occurs at the temperature much lower than the 
temperature of the magnetic transition. Normally differ- 
ent elements of the crystal structure are responsible for 
ferromagnetism and superconductivity, and spontaneous 
magnetization (ferromagnetic order parameter) is related 
to spin. Later we call them spin superconducting ferro- 
magnets (spin SFs). (ii) The magnetic and the SC transi- 
tions occur simultaneously. This can take place in uncon- 
ventional superconductors with triplet Cooper pairing. 
An example of them is strontium ruthenate Sr2Ru04 
[1, 0, Q- The theory connects spontaneous magnetiza- 
tion in this material not with spin but with the orbital 
intrinsic angular moment of the p-wave Cooper pair with 
the wave function in the momentum space proportional 
to p x + ip y (chiral p-wave superconductivity). We shall 
call them orbital superconducting fcrromagnets (orbital 
SFs). 

Whereas proof of superconductivity in SF materials 
is quite straightforward, a clear-cut detection of the fer- 
romagnetic order parameter is much more problematic. 
The internal magnetic field is screened out by the SC 
Meissner currents and can be present only near sample 
borders and defects, in particular, domain walls (DWs). 
This strongly suppresses the stray magnetic fields around 
the sample, which are most convincing evidence of ferro- 
magnetism. Especially worrying is situation with stron- 
tium ruthenate S^RuO/i, where Kirtley et al. [8| could 
not detect any stray field from DWs or sample edges at 
all. This is a challenge for the theory and for the very sce- 
nario of chiral p-wave pairing. Difficulties with direct de- 
tection of ferromagnetism coexisting with superconduc- 
tivity lead to the question whether one may use the term 
ferromagnetism at all. Indeed in the literature on uncon- 
ventional superconductors sometimes they prefer to tell 
about superconductivity with broken Time-Reversal Sym- 



metry (TRS). We used this name in the title of the paper 
following this more cautious semantics though one can- 
not imagine broken TRS without at least some features 
of ferromagnetism (see further discussion). 

Among possible explanations why experimentalists 
cannot see stray fields from DWs is the presence of do- 
main structure with a period essentially smaller than 
a distance between a sample surface and a probe used 
by experimentalists. There were some experimental ev- 
idences of domains in SFs both in the spin [9( and the 
orbital SF [lfj . The theoretical investigations of the do- 
main structure in SFs were restricted with the case of 
zero external magnetic field [ll], [ID, [TH, [l4| • One must 
discern two possible types of equilibrium domain struc- 
ture. The first one is well known for normal ferromag- 
nets [lol ]. The domain structure results from competi- 
tion between the energy of DWs and the magnetostatic 
energy of stray fields generated by the magnetic flux exit- 
ing from the sample surface. The period of the structure 
depends on the shape and the size of the sample going 
to infinity when the sample size grows. One can call 
these domains extrinsic ferromagnetic domains. Since in 
SFs the Meissner effect expels the magnetic field, it is 
impossible to benefit from decreasing the bulk magne- 
tostatic energy in comparison with the DW energy, and 
extrinsic domains cannot appear at equilibrium [121 ] . But 
also long ago there was known another type of domains, 
which decrease the bulk magnetostatic energy at the ex- 
pense of destroying the Meissner state [ill EE EH ■ The 
size of these domains is roughly of the order of the Lon- 
don penetration depth A and does not depend on either 
shape or size of the sample. Strictly speaking the state 
with this domain structure at the macroscopic scales is 
not ferromagnetic but antiferromagnetic though with a 
rather large period. We shall call such a state cryptofer- 
romagnetic, the term introduced by Anderson and Suhl 
16] for another model of ferromagnetism and supercon- 
ductivity, in which crystal anisotropy was neglected and 
spiral structure appeared instead of domains. 

In present publication we extend the theory of intrinsic 
domain structure (cryptoferromagnetic state) on nonzero 
external magnetic field and analyze competition of the 
cryptoferromagnetic state with the pure Meissner state 
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and the mixed state with vortices. This allows to ob- 
tain the full phase diagram of both spin and orbital SFs. 
We demonstrate that the measurement of magnetization 
curves in various areas of the phase diagram can provide 
evidence of ferromagnetic or cryptoferromagnetic order 
in superconductors with broken TRS. 

Let us consider a stripe domain structure with 180° 
DWs in a sample subjected to external magnetic field 
H = (0, H o ,0). The DWs are parallel to the yz-plane 
separating domains with alternating magnetization M = 
(0,±M ,0) along the +y or — y direction. Since the H 
orientation is preferable the width d-\ of domains with 
the magnetization M parallel to Ho (t-domains) exceeds 
the width d[ of the domains with M antiparallcl to Ho 
(|-domains). We restrict ourselves to the simplest case 
when the London penetration length A exceeds the DW 
thickness. Then the surface energy a and the internal 
structure of DW is not affected by fields and currents at 
scales of A. 

The Gibbs potential inside domains is 
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where h is the magnetic field, and the electric current j 
is connected with the magnetic field h via the Maxwell 
equation V x h = (4n/c)j. Variation of the Gibbs po- 
tential yields the magnetic field hj^ = (0,^^,0) in the 
t-domains and |-domains: 



h u = {H ± 4ttM ) 



cosh (x/X-^ d ) 
cosh£ u 



(2) 



where x is the distance from the DW and = dj.j/2A 
are reduced domain widths. 

Application of the Gibbs potential Eq. (Q} to orbital 
SFs requires some comments. As shown in Ref. [l7j . for 
orbital ferromagnetism related to the intrinsic angular 
momentum of Cooper pairs the spontaneous magnetiza- 
tion cannot be defined unambiguously and therefore the 
Landau-Lifshitz theory of ferromagnetism [l5[ based on 
this definition is not valid. Nevertheless, interaction of 
magnetization currents inside narrow DW with the mag- 
netic field can be reduced to the expression looking like 
the standard Zeeman energy — h • Mo- However, here 
Mo is not a magnetic moment inside the domain but is 
defined so that 8ttMo would be the jump 8ttMq of the 
magnetic field on the DW [see Eq. So "magnetiza- 
tion" Mo is determined by the DW structure and cannot 
be used for other phenomena connected with ferromag- 
netic ordering, e.g., analyzing the magnon spectrum [17| . 

Substituting Eq. @ into Eq. |T]), adding the sur- 
face energy a of DWs, and averaging over the domain- 
structure period d = d-\ + we arrive to the following 
expression for reduced energy density £ — Q/2ttMqV (V 
is the sample volume): 

_ 2w- (1 + /i ) 2 tanh£ T - (1 - fe ) 2 tanh^ 



Here ho = Ho/^nMo and w = a/AnM^X are dimension- 
less parameters. If h = and £j = £j Eq. © coincides 
with the free energy density of Krey [1JJ . 

Minimization of energy density Eq. ([3]) with respect to 
£f and £j yields the system of two nonlinear equations 
for £ T and ^ : 
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The magnetic induction B = (h) is determined by re- 
duced magnetic induction b = B/AirMo'. 



1 d£ _ (1 + h ) tanh£ T - (1 - ho) tanh^i 

2dho~ e T + ei 
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Fig.[lja) shows the phase diagram in the plane w 2 ~ho- 
The area of the cryptoferromagnetic state is restricted by 
two lines where the phase transition between the crypto- 
ferromagnetic and the Meissner states occurs: 

1. The line 6 = 0, which corresponds to the limit 
— > oo. The values of w and ho on this line are con- 
nected by the relation w = 1 + h^ . 

2. The line on which domains with magnetization op- 
posite to the external magnetic field vanish, £j = 0. The 
equation describing this line is 



U'c 



2 " 1 - Vh c 
the critical size of the "["-domain being 

l + xfhr 



In- 



1- Vh c 

The magnetic induction on the critical line is 

1 + sfhe. 



b r = 2\/h r \ In 



1-Vhf 



(6) 
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Let us consider the left lower corner of this diagram 
where w <Cl and ho <C 1. The critical parameters on the 
line £| = as functions of w are 



J(3H 2/3 , 



Aside from the critical line Eqs. (J3]) yield: 

ho = 



(9) 
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For small ho <C h c one can solve Eqs. (JTUJ) analytically: 

ho 
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Let us consider the magnetization curve in the crypto- 
ferromagnetic state. The linear magnetic permeability is 
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FIG. 1: (color online) Phase diagram for various values of the 
reduced lower critical field h&: (a) h c i — > oo; (b) h c \ = 2; (c) 
h c i — 0.83. The lighter (yellow) shaded area is the cryptofer- 
romagnetic state. The darker (blue) shaded area is the mixed 
state. The rest is the Meissner state. The horizontal and the 
vertical arrows show the processes of cool-down across the SC 
critical temperature of spin and orbital SFs respectively. 



determined from two relations connecting \i and w with 
the period £ = £ T + £ x w 2£ T 

db coth£ , ^ £ , . 

^=^ = V' W = t ™ H -^ (12) 

In the limit w — > the magnetic permeability is diver- 
gent: jit w (2/3w) 2 / 3 . The whole magnetization curves 
b(ho), which were calculated numerically, are shown in 

Fig. EDO. 

Fig. I2jb) shows the dependence of the period £ = £j + 
and the difference of the two domain widths S = — 
as functions of the reduced external magnetic field ho- 

Up to now we ignored the possibility of the transi- 
tion to the mixed state assuming that the first critical 
magnetic field H cl = <f> lnn/AirA 2 essentially exceeds 
the characteristic fields of the cryptoferromagnetic state. 
Here k = A/£o is the ratio of A to the coherence length 
£o- Let us now take into account this possibility. Both 
the field H c \ and the parameter w depend on the pene- 
tration depth A, and it is useful to introduce the reduced 
first critical field h cl = H c i/AttMqw 2 — $qMq Inn/a 2 , 
which does not depend on A. The reduced free energy of 
the mixed state with respect to the energy of the Meiss- 
ner state is 




where the field H* inside the mixed state differs from H c \ 
by another logarithm factor, but we neglect it assuming 
H* H c i. The phase transition to the mixed state may 
occur either from the Meissner state being determined by 
the condition £ m = 0, or from the cryptoferromagnetic 
state crossing the critical line on which £ m = £ . At zero 
external field h$ and small w the phase transition be- 
tween the mixed state and the cryptoferromagnetic state 
occurs at w m ~ \/3/(2h c i) 3 ^ 4 . Thus whatever large h c i 
is, in the left lower corner of the phase diagram there is 
always the spontaneous vortex phase, i.e., the mixed state 
without external magnetic field. The full phase diagrams 
at two finite values h c i —2 and 0.83 are shown in Figs, 
[ljb) and (c). The cryptoferromagnetic state disappears 
from the phase diagram at h c i < 0.5. 

Now let us analyze the phase transformations in the 
process of cooling down below the SC critical tempera- 
ture. This process is different for spin and orbital SFs. 
In the case of spin SFs, when the magnetic transition 
occurs at much higher temperature, one may neglect 
temperature dependence of Mo, a, and h c i. Then only 
w 2 oc 1/A 2 oc t depends on relative temperature differ- 
ence t = (T c — T)/T c . On the phase diagrams of Figs. 
[TJb) and (c) the state moves along straight lines paral- 
lel to the horizontal axis w 2 . From these figures it is 
evident that just below the critical temperature the sys- 
tem enters the mixed state. At further cooling down the 
system crosses to the Meissner state either directly or 
through the area of the cryptoferromagnetic state. For 
orbital SFs the cooling process occurs differently. In 
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FIG. 2: (color online) (a) Magnetization curves and (b) mag- 
netic field dependencies of f = £(/io) (solid line) and S = 8(ho) 
(dash line) are shown at different values of parameter w. Ver- 
tical lines correspond to critical fields above which the intrin- 
sic domain structure collapses. 



this case the "magnetization" Mq ~ $/A 2 oc r, and the 



DW surface energy is a product of the condensation en- 
ergy H%(t) ~ [$o/-M 7 ")£o(' 7 ")] 2 and the coherence length 
Co(t) ~ fo/V?: o- - T 3 / 2 $2/A2,e . Here A and £ are the 
penetration depth and the coherence length at zero tem- 
perature. Then the parameters w 2 ~ h~i ~ (Ao/£o) 2 do 
not depend on temperature whereas the reduced mag- 
netic field does: ho — Hq/AttMq cx 1/t. Thus in the 
field-cooling process the state moves along vertical lines 
on the phase diagrams in Figs. QJb) an d ( c )- However, 
as pointed out above, the cryptoferromagnetic state can 
compete with the mixed state only if h c \ is high enough. 
Since h c \ ~ (<^o/Ao) 2 , this requires Ao not large com- 
pared to £o- According to 0] the ratio of Ao = 190 nm 
t° Co = 66 nm is not too high indeed. But this means 
that the DW thickness is not so small compared to A as 
assumed in our analysis. Therefore for orbital SFs our 
analysis can provide only a qualitative but still credible 
picture of the phase transformations. 

In conclusion, we analyzed the conditions for appear- 
ance of intrinsic domains in superconductors with bro- 
ken time reversal symmetry, in which superconductivity 
coexists either with spin ferromagnetism or with ferro- 
magnetism originated from the intrinsic angular momen- 
tum of Cooper pairs (chiral p-wave superconductors like 
Si'2Ru04). Since these domains strictly speaking corre- 
spond not to a ferromagnetic but a globally antiferro- 
magnetic state, the state was called cryptoferromagnetic. 
We considered competition of this state with the Meiss- 
ner and the mixed states and found the phase diagram. 
This phase diagram can be checked with detailed mea- 
surements of magnetization curves. 
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